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Abstract 

We discuss a 3-flavour lattice QCD action with clover improvement 
in which the fermion matrix has single level stout smearing for the hop- 
ping terms together with unsmeared links for the clover term. With the 
(tree-level) Symanzik improved gluon action this constitutes the Stout 
Link Non-perturbative Clover or SLiNC action. To cancel 0{a) terms 
the clover term coefficient has to be tuned. We present here results of 
a non-perturbative determination of this coefficient using the Schrodinger 
functional and as a by-product a determination of the critical hopping 
parameter. Comparisons of the results are made with lowest order pertur- 
bation theory. 
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1 Introduction and 0(a) improvement 



When constructing a lattice QCD action, even the simplest gluon action has only 
O(a^) corrections. The naive quark action also has O(a^) corrections, but suffers 
from the 'doubling problem' describing 16 flavours in the continuum limit. A 
cure is to add the Wilson mass term, so 15 flavours decouple in the continuum 
limit, but the price is that there are now 0(a) corrections (and also loss of chiral 
invariance), so that for example for a ratio of hadron masses 

I^=ro + ari + 0(a2). (1) 

The Symanzik approach is a systematic improvement to 0(a") (where in prac- 
tice n = 2 for the fermion action) by adding a basis (an asymptotic series) of 
irrelevant operators and tuning their coefficients to remove completely 0(a"~^) 
effects. Restricting improvement to on-shell quantities the equations of motion 
reduce the set of operators in both the action and in matrix elements. Indeed, 
for 0(a) improvement of the fermion action only one additional flavour-singlet 
operator is required 

Cdover acsw ^ q{x)a ^^F^y{x)q{x) , (2) 

the so-called 'Sheikholeslami-Wohlert' or 'clover' term, [I]. So if we can improve 
one on-shell quantity this then fixes function of the lattice spacing a 

or equivalently of the bare coupling g'^, so that all other on-shell quantities are 
automatically improved to 0(a), i.e., we now have 

^ = ro + 0(a^). (3) 
rriH' 

A non-perturbative determination of Csw will be the main goal of this paper, the 
general approach being described below. 

Matrix elements still require additional 0(a) operators, for example for the 
axial current and pseudoscalar density, [2fi 

A^, = {l + bAamg){A^ + CAad;fP) 
V = (l + 6pam,)P, (4) 

(for mass degenerate quarks) with 

^/^ = llt^l^q , P = qim , (5) 



^We implicitly distinguish between quark flavours in operators, i.e. consider non-singlet 
operators. 
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which require additional 6a , ca and hp improvement coefficients. An easily de- 
termined quantity is the quark mass computed from the PCAC WI relatioio, 

~ 2(P(xo)0) • 

Choosing different operators, O, gives different determinations of the quark mass 
m^^^*\ i = 1, 2 with different lattice artifacts. If the quark mass is improved 
then its errors are O(a^). So we can determine the 'optimal' Cgw improvement 
coefficient by tuning until 

(This is equivalent to considering the renormalised quark mass 

_ Za{1 + bAarrig) 



Zp(l + bpam, 



qj 



In general the Ba, bp coefficients do affect considerations of 0(a)-improvement. 
However, here one imposes a condition at fixed bare parameters {gQ,mg) which 
means that the factors drop out.) Practically, how this is achieved will be dis- 
cussed in this paper after the action is introduced. 

This paper is organised as follows. In the next section, section [2|, the action is 
given and in the following section the Schrodinger functional is briefly discussed, 
together with the general procedure for determining the optimal Cgw and optimal 
critical hopping parameter, k^. Section H] gives some lattice details for a series 
of simulations at various coupling constants, which after suitable interpolations 
leads to this determination. Section [5] then discusses possible finite size effects in 
the results. Results are collected together in section [6] and a polynomial interpo- 
lation (in the coupling constant) for both Csw and are given, together with a 
comparison with the lowest order perturbation result. Finally in section [7] some 
brief conclusions are discussed. Tables of the raw results are given in appendix Rl 



2 The SLiNC action 

We shall consider here nf = 3 flavour stout link clover fermions - SLiNC fermions 
(Stout Link Non-perturbative Clover). In a little more detail, for each flavour, 

Sf = 

H \ - ^)U^{x)q{x + afi) - g(x)(7^ + l)t/^(x - afi)q{x - afi)] 

+q{x)q{x) - i/^ac^^„^g(x)o-^^F^l,(x)g(2;) \ . (9) 

/lU ) 

^dj^ is the symmetric lattice derivative, f){x) = [/{x + afi} — f{x — ajl)]/(2a), and (no 
fl summation), {df,"^ f){x) = [f{x + afi) - 2f{x) + f{x - afl)]/a^ = {d^d^f){x) + 0{a^). 
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Rescaling the quark fields q ^ q/ v2k gives the quark mass where 

1 1 \ 




(10) 



which is proportional to the PCAC quark mass, m^' . The loss of chiral invariance 
means that for a given Csw a critical hopping parameter, Kc(cs^) has now also to 
be determined. 

The hopping terms (Dirac kinetic term and Wilson mass term, i.e. those terms 
involving a k) in eq. ([9]) use a once stout smeared link or 'fat link', [3], 

Uf, = exp{iQf,{x)}U^{x) 

Q^{x) = ^ [VU^ - UV^ - iTr(\/f/t - UV^)] , (11) 

(V^ is the sum of all staples around f/^) while the clover term remains built from 
'thin' links - they are already of length 4a and we want to avoid the fermion 
matrix becoming too extended. Smearing is thought to help at present lattice 
spacings by smoothing out fluctuations in the gauge fields slightly and so reducing 
the condition number and also to avoid a near first order phase transition. The 
critical kappa in eq. (|T0|) corresponds to an additive mass renormalisation 



mc(c,^) = — — r - • (12) 

It is known that with a combination of link fattening and increase of the clover 
coefficient, it is possible to reduce this mass term [H El E]- The stout variation 
is also analytic which means that the derivative in the gauge group can be taken 
(so the force in the Hybrid Monte Carlo, or HMC, simulation is well defined) and 
perturbative expansions are also possible, [7]. 

To complete the action we also use the Symanzik tree-level gluon action 

Sg = (13) 

Y ^ReTr(l - f/p,,^„,,,J + Cl Y ^^eTT{l-Uaectangle) \ , 

Plaquette Rectangle ) 




together with 

20 ^ A R 6co 10 

Cn = — , Cl = and p = ^5- = ^ . (14) 

12' ' 12 ^ gl gl ^ ' 

While this gluon action has elements of higher order improvement, namely O(a^), 
this is not the reason that it is used here. (The best we can hope for the fermion 
action is O(a^) improvement.) Again we wish to move the action away from 
a nearby first-order phase transition occuring when using the standard Wilson 
action (i.e. cq — 1, ci — 0), [8] by using a slightly extended action. Different 
values of cq and ci can be and have been used in the literature to address this 
problem, e.g. [8]. 
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3 The Schrodinger functional 



The ALPHA Collaboration determined the improvement coefficients by means of 
the 'Schrodinger functional', [9l [TOl [TTl [2] . Some numerical results for c^^ for the 
quenched case {nj = 0) were given in [121 [13]) fo'^ '^Z ~ ^ flavours in [H] and for 
Uf = 3 flavours in [151 [16l [17] . In this approach Dirichlet boundary conditions are 
applied on the time boundaries to the fields. For the gluon fields, fixing them on 
the boundary is then equivalent to inducing some classical background field about 
which they fluctuate. It is simplest to consider spatially constant colour diagonal 
fields, corresponding to a constant chromo-electric background field. Concretely, 
we consider a x T lattice (with T = 2L) and take the background field to be 



Uq{x,xo) 



exp(-^|[xoC(2) + (T 



Xq 



1)' 



with 



\ 


V 0« 







and 



4". 



i>r,0, 



in 



(15) 



(16) 



;i7) 



and fix the boundary values a posteriori. As we have an extended gauge action 
(rather than the simpler Wilson gluon action), we fix two values at each double 
boundary layer and so we choose, following [isj^ . from eq. (ITS]) at Xq = —a, 
(lower boundary) and similarly [/^ at a;o = T — a and T (upper boundary). 
The 'bulk' of the lattice is thus from xq = to xq = T — a. Additionally the 
weight factors for the gluon loops in eq. ([14]) must be appropriately chosen on the 
boundary for 0(a)-improvement. Classically these weight factors are not difficult 
to find, however a full non-perturbative determination would be difficult. But 
away from the boundaries, they only affect the local PCAC relation to O(a^) and 
so are not essential for the determination of the optimal Csw, and so it is sufficient 
to use the classical values. 

The fixed boundary quark fields, p, p (taken as zero here) make simulations 
with rriq ~ with no zero mode problems possible. They are specified on the 
lower inner boundary and upper inner boundary from 



Po+g(^,o) 

g(x,0)Po" 
PoQ{x,T- a) 
q{x,T-a)P(^ 



p«l 
p(^)( 
p(^)| 



X) 



x) 



x) 



x) 



^An alternative procedure using single layer boundaries is given in |19j . 
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where Pq is the projection operator defined by 

Po^ = i (1 ± 7o) . (19) 

These projections are necessary for consistency, p, p can be taken as sinks and 
sources respectively to build operators for correlation functions. For example 
here we can take at the lower inner boundary xq = [i = 1) and upper inner 
boundary xq = T — a {i = 2) the operators 



So we can investigate PCAC behaviour at different distances from the boundaries. 
In a little more detail, following eq. ([6]), we first set 



^W^^^) = Z ; " > ^'"Kxo) = a \ Jl ' , (21) 
where 

(^o) = -^^(Ao(xo)0«) 
n)-\ 

ff{T-x^) = +^l-(Ao(xo)0(2)), (22) 



nj-l 



and 



f?ixo) = -^l-(P(a;o)0«) 



n}-l 



ff{T-xo) = -^l-(P(a;o)0(2)). (23) 



Then redefine the quark mass slightly, but which coincides to O(a^) for the im- 
proved theory 

M»(xo,yo) = r»(a:o) + CA{yo)s^Hxo) , Uvo) = -^TTT^, ^ , (24) 



which eliminates the unknown ca in the determination of the quark mass, [12] 
and replaces it by an estimator, c^. Improvement is defined when 

(M, AM) = (0,0) , (25) 

where 

M = M(^) , AM = M(^) - M(2) , (26) 
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are chosen at some suitable xq, [12]. This gives the required optimal Csw and Kc, 
which we will denote by a star: c*^ and k*. Conventionally, we choose 

M = M(^) (T/2, T/4) , AM = M(i) (3T/4, T/4) - M^^) (3T/4, T/4) . (27) 

There are small changes due to the finite volume used, so eq. (l25ll becomes 

(M, AM) = (0, AM*"*^") , (28) 

where AM*'"'^'^ is the tree-level (i.e. = 0, Csw = c^sw'^ = 1) value of AM|m=o on 
the X T lattice. This ensures that Csw 1 exactly as (3 —>■ oo. For a = 0, 
the analytic result on a A^^^ x 2Ns = 8^ x 16 lattice (where L = aN^) is 0.000277, 
|12] . Carrying out the interpolation procedures outlined in the next section for 
a free configuration, with background field given by eq. ( fTSl) yields 0.000271. 
For the stout smearing used here (see next section) we find this is reduced to 
^j}^tree _ 0.000066 and so we have neglected AM*'''^*^ in the following and simply 
used eq. (!25|) . 



4 The lattice simulation 



The 3-fiavour lattice simulation used the Chroma software library, [20], the 
Schrodinger functional details following [18]. Results were mostly generated on 
Ng X 2A^s = 8^ X 16 lattices, together with some additional 12^ x 24 lattices, using 
the HMC algorithm together with the RHMC variation, [21], for the 1-fiavour. 
A mild smearing of a = 0.1 was used. A series of simulations were performed 
(typically generating 0(3000) trajectories for the 8'^ x 16 lattices and 0(2000) 
trajectories for the 12'^ x 24 lattices), quadratic and then linear interpolations of 
the (M, AM) results being used to locate the optimal point, (0, 0) as described 
below. Some further details and tables of the results are given in appendix |Al 
(Preliminary results were given in [22].) 



We have a two-parameter interpolation in Csw and k which is split here into two 
separate interpolations. First plotting AM against M and then interpolating to 
M = for fixed Csw gives a critical n namely Kc{csw), 

AM(c,^, k)\m=o = AM(c,^, Kc(cs^))|m=o = AM(c,^) . (29) 

In Figs. [H [2] we plot AM versus M for various Cgw values for the 8^ x 16 lattices 
and in Fig. [3] the results for the 12^ x 24 lattices. 

These graphs are the fundamental plots requiring high statistics as AM is the 
difference between two different Ms. As there are always 4 (or more) points for 
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Figure 1: AM against M for /? = 5.10, 5.25 (upper left, right pictures respectively) for 
P = 5.50, (3 = 6.00 (middle left, right pictures respectively), and for /3 = 6.50, /3 = 7.20 
(lower left, right pictures respectively), together with quadratic interpolations to M = 
(the open symbols). 
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Figure 2: AM against M for /? = 8.00, 10.0 (upper left, right pictures respectively) 
and for /? = 14.0, (lower picture), together with quadratic interpolations to Af = (the 
open symbols). 
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Figure 3: AM against M for (3 = 5.50, 6.00 (left, right pictures respectively) on a 
12^ X 24 lattice together with quadratic interpolations to M = (the open symbols). 
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Figure 4: AM at M = against Cgw for various values of /? (filled circles) together 
with linear interpolations to AM = (open circles). The left plot shows the 8^ x 16 
results while the right plot shows the 12^ x 24 results. 



each graph a quadratic fit is made and the value of AM is determined where M 
vanishes. 

These values of AM(cs^) for each (3 value are then plotted against Csw as 
shown in Fig. H] together with linear fits. The point where AM{csw) vanishes 
gives c*^. This gives values of 
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We postpone a discussion of possible finite size effects until section [51 

From Fig. HJ we see that linear fits even for four points (the (3 = 7.20, 6.00, 
5.25 results) show very little curvature, so that we may write, [H] 

AM(c,^) = u (c,^ - clJ , (31) 

with the gradient, u, a slowly varying function of qq. To test this we note that 

5AM(c,^ 



(32) 
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so a fit to the gradients in Fig. H] (for tlie 8^ x 16 lattices) yields an estimate for 
uj. We find that uo is constant with an approximate value of —0.018, although for 
the largest values of there are deviations from this. 



4.2 



til 



A similar procedure yields k*: plotting M against 1/k and interpolating quadrat- 
ically to M = for fixed Csw gives the critical k, denoted by Kc(cs«,). Then sub- 
sequently plotting iS.M{csw) against 1/ Kc{csw) and interpolating using a linear fit 
to AM = gives k*. 

We first plot M against 1/ k for the 8^ x 16 results in Figs. El Eland for the 
12^ X 24 results in Fig. El 

Note that to produce these graphs should not require high statistics as it does 
not involve AM. (Although these are not the fundamental graphs they are also 
useful in helping to determine the various {csw, i^) values for the runs.) 

These ^M{kc) are then plotted in Fig. [8] again with a linear fit. Where AM 
vanishes gives k*. For legibility the results have been split into sub-graphs. We 
see that is a non-monotonic function of (3. 

We find results of 
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12^ X 24 



As a consistency check the alternative plot of Csw against 1/kc is shown in 
Fig. where Cgw is plotted against 1/kc{csw), again with a linear fit between 
the points. The optimal values of Cs^, namely c*^, taken from the previous fits 
as given in eq. (l30l) are shown as dashed horizontal lines, the intersection with 
the 1/kc curves then giving the optimal critical values of Kc, namely /t*. These 
are denoted in the figure as open points. As a comparison, the results from 
the previous determination of k*, eq. (155]) . are also shown as vertical lines. We 
see good agreement between the different determinations of k*, which indicates 
that the fit procedure adopted here gives consistent results for both c*^ and k*. 
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Figure 5: M against 1/k for /3 = 5.10, 5.25 (upper left, right pictures respectively) for 
(3 = 5.50, P = 6.00 (middle left, right pictures respectively), and for (3 = 6.50, f3 = 7.20 
(lower left, right pictures respectively), together with quadratic interpolations to M = 
(the open symbols). 
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Figure 6: M against 1/k for /3 = 8.00, 10.0 (upper left, right pictures respectively) 
and for /? = 14.0, (lower picture), together with quadratic interpolations to M = (the 
open symbols). 
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Figure 7: M against 1/k for /? = 5.50, 6.00 (left, right pictures respectively) on a 
12^ X 24 lattice together with quadratic interpolations to M = (the open symbols). 



13 



• p=14.0 



0.000 ^> 







• |3=5.50 

• §=6.00 




• p=6.50 T 

I -"i 




\ 1 ' 


\ 


-f T 






\ 


\ i 


1* 



i. 



3.0 8.2 6.4 8.6 

1/K. 



3.0 8.2 8.4 6.6 

1/k 



Figure 8: Results of AM{kc{csw)) versus 1/kc together with linear fits. The open 
circles give the optimal critical KcS, i-e. the k*s. The two left plots show the 8^ x 16 
results while the two right plots compare the /5 = 5.50, 6.00 8^ x 16 results with the 
12^x24 
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Figure 9: Results of (filled circles) versus 1/kc together with linear fits. The 



optimal Csw, c*^, from eq. (j30j) are shown as dashed horizontal lines. The open circles 
are the intersection of the linear fits with these horizontal lines and give an alternative 
determination of the optimal critical Kc, k*, which are to be compared with the results 
of eq. (j33p shown as vertical lines. 



14 



Finally note that plotting the Uf = 2 flavour results would yield a similar curve 
to Fig. El 

For future reference (in section 16.21) as the fits in Fig. [9] are all linear then we 
write 

- = — + dicsw - clJ , (34) 
with a measured coefficient d{go), 
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(35) 



5 Finite size effects 



There are (small) ambiguities due to the finite volume used. In an infinite vol- 
ume we expect ©(aAqci,) contributions (in the chiral limit, otherwise there are 
also extra 0{amg) terms) due to the different boundary conditions or operators 
chosen. In a finite volume there are additional 0{a/L) terms. Thus might expect 
asymptotically, following |T6], 



cL(^o, L/a) = cl^igo, oo) + cl- + ca aAocD + ... ■ 



(36) 



The terms proportional to aKqcD vanish as a 



(or gl) 



and represent the 



ambiguities in the different definitions of M. For a physical quantity Q, then 



Q = Qia) + qLicUgo,L/a)-c*Jgo,oo))aAQao + 0{a'') 



SW \ 
2\ 



= Q(a) + QlCl j aAgcn + 0{a') . 
The correction term may be re-written as (where L = aNg 



(37) 



qlCl — aAoco = -j^ a A 



QCD 



(3^ 



Potentially this might mean that Q is no longer 0{a) improved for simulations 
where c*^ has been determined on a fixed lattice size, Ng. However it is likely 
that the unknown coefficients qi and cl are small and coupled with the Ns factor 
in the denominator, this is then expected to be a small effect. 
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Figure 10: c*^ against l/Ng (left picture) and k* against l/Ng (right picture) for 
(3 = 5.50, 6.00, filled circles. Also shown are constant fits (dashed lines) together with 
the extrapolated values (open circles). 



To avoid this altogether we can either keep L fixed in physical units as a — 
(the 'constant physics condition') so 0{a/L) — > 0, or alternatively simulate for 
several values of Ns and extrapolate to iV^ oo. The 'Poor man's solution' is 
to evaluate at large (3 ^ oo (i.e. on a free configuration for iV^ = 8 here) and 
subtract this result. Practically, following the same procedure as in section 14.11 
we have found that for Csw this 0{1/Ns) term (for A^^ = 8) is negligible. 

As noted previously we have also performed additional simulations on larger 
lattices 12^ x 24 for (3 = 6.00, 5.50 in order to discuss finite lattice size corrections. 
The results are plotted in Figs. [H [7] and compared with the 8^ x 16 results in 
Figs, m [HI At tree level we have, [23], 

AM'^'' = k {c';:' - l)y + . . . , (39) 

which would indicate that for larger Ng then AM becomes smaller, with the 
consequent noise/signal ratio becoming worse. Indeed this is seen in our results, 
with the 12^ x 24 data being more bunched together in Fig. [3] than for the 
corresponding 8^x 16 data in Fig.[TJ This may be mitigated somewhat by choosing 
a larger range of Csw due to the linear nature of the data as seen in Fig. H] and 
eq. (I3T1) . For P = 6.00 we have increased the number of CswS used in the analysis. 

In Fig. [TU]we plot c*^ and k* against l/Ng. For both /? = 6.00 and 5.50 there 
seems to be small finite size effects for c*^. For k* this is also the case for (3 = 6.00, 
while for j3 = 5.50 the situation is perhaps a little less clear-cut. However there 
is no systematic trend in the data and a constant fit always lies within the error 
bars of the data. So although we cannot come to a definite conclusion, there do 
not seem to be large finite volume effects, i.e. cl appears to be small in eq. fl5B]) . 
So in eq. fl?r|) we only expect small violations of 0{a) improvement. We shall, in 
future, just consider the 8^ x 16 data. 
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6 Results for c^^ and 

6.1 Perturbative results for c^^ and 

Before giving the non-perturbative results for c*^ and k* we first recapitulate the 
perturbative results. The lowest order perturbative limit has been computed for 
both c*^ and n*, f7\- For c*^ we have 

cl^{go) = 1 + (0.196244 + 1.151888a - 4.2391365a^)^o , (40) 

where a is the stout smearing parameter, set equal to 0.1 here. This gives 

CsJgo) = 1 + cigl , ci = 0.269041 , (41) 

i.e. the smearing parameter has increased the value of c*^ (for a = 0, we have 
ci = 0.196244). For Kc{csw,go) we have 

Kc(c,^, go) = + (0.0853699 - 0.961525a + 3.55806a^ (42) 

8 '- ^ 

-(0.025221 - 0.0787379a)c,^ - 0.00984224cy gl 
giving for a = 0.1 

t^cicsw, go) = lh + (0.024798 - 0.0173472c,^ - 0.00984224cL) ^ol > (43) 



and finally for Cgw = cfw"^ = 1, 

<{go) = ^ fl + hgl] , h = -0.002391 . (44) 

(Note that the result for k,c{csw, go) is more general than the one given in [7] when 
only the result for Csw = 1 was given.) 



6.2 Non-perturbative results for c*^, and 

The results for c*^ and against are plotted in Figs. [TT1[T2] respectively in the 
range (3 > 5.10. The lowest order perturbative limits are also shown, eqs. (HTi) 
and dSD. 

An interpolation between the numerically determined points is also shown. 
For both c*^ and k* a 5th order polynomial in g^ proved sufficient. (These 
interpolation functions are constrained to reproduce the perturbative results, in 
the P CO limit and therefore, they have four free parameters.) For c*^{go) we 
write 

cLigo) = 1 + Cl^o + C2^o + CS^O + C4^o + ^5^0° , (45) 
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Figure 11: c*^ against Qq for various vaiues of /? (circles), together with a polynomial 
interpolation (line). Also shown is the perturbative result. 



and find 



C2 

C4 
C5 



+0.29910 
-0.11491 
-0.20003 
+0.15359 



while for K^go) we write 



(46) 



10 



(47) 



and find 



k2 

h 

k5 



+0.0122470 
-0.0525676 
+0.0668197 
-0.0242800 



(48) 
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Figure 12: k* against gQ for various values of f3 (circles), together with a polynomial 
interpolation (line). Also shown is the perturbative result. 



These give for the specific (3 values used here 



3.306 


[3 = 


5.10 


' 


116262 


f3 = 


5.10 


3.021 


P = 


5.25 





118424 


P = 


5.25 


2.653 


P = 


5.50 





120996 


P = 


5.50 


2.179 


P = 


6.00 





123751 


P = 


6.00 


1.907 


P = 


6.50 < = < 





124870 


P = 


6.50 


1.692 


P = 


7.20 





125328 


P = 


7.20 


1.560 


P = 


8.00 





125314 


P = 


8.00 


1.407 


P = 


10.0 





124979 


P = 


10.0 


1.279 


P = 


14.0 


. 


124783 


P = 


14.0 



which are to be compared with the numerically determined values. The errors 
for c*^ from the fit are estimated to be about 0.4% while for k* we have 0.02% 
at /3 = 14.0 rising to 0.15% at /? = 5.10. 

These smooth fits between the points give estimates for c*^ (and k*) which 
could be used in the action for future generation of configurations. 

For c*^ the polynomial only tracks the perturbative solution for small values 
of (7q. This is perhaps not surprising as the tadpole improved, TJ, estimate is 

= Uq^^ /uq, [7J, which is to be compared with the unsmeared case of cj^ = 1/mq 
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(S) 

where uq is the average plaquette value and Uq is the smeared value. As smearing 
increases the plaquette value this indicates that c*^ can be large. For on the 
other hand as = l/(8'Uo ) we expect that it is ~ 1/8. This is true for 
reasonably fine lattices, however k* does begin to decrease for larger values of Qq. 
For Uf = 2 the same phenomenon occurs: for larger Qq, k* begins to decrease 
(after initially increasing). 

As a further consistency check on the results, we can investigate the gradient 

'-'SW \( 



From eq. (IMI) we have 
d{l/K, 



d . 



(50) 



as the fits in Fig. [9] are linear, where d is given in eq. (l35l) . Perturbatively we 
have from eq. (H3|) . 



8 [0.037032 + 0.019684(c,^ - 1)] . 



(51) 



As qq increases Csw increases, so not only do more terms in this expansion become 
important, but the coefficient of the leading term increases as well. For Csw = 



^tree 



1 we have the leading order perturbative result. 



d^gl , di = 0.296253 . 



(52) 



In Fig. [13] we plot d{l/Kc)/da 
polynomial in Qq, 



against Qq, together with a 5th order 
diQo + d2gl + d^gl + d^gl + d^gl^ , (53) 



and find 



d2 


+0.4180 


ds 


-0.7232 


di 


+0.4739 


d5 


-0.0919 



(54) 



The results follow a smooth curve. 



7 Conclusions and Discussion 

Non-perturbative 0{a) improvement is a viable procedure for (stout) smeared 
actions with typical clover results being obtained. (Other recent results for 3 
fiavours are given in [151 [13 [13 •) Using the Schrodinger functional method 
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Figure 13: d{l/ Hc) / dcsw\ci^ against for various values of /3 (circles), together with 
a polynomial interpolation (line). Also shown is the perturbative result. 



we have determined the optimal clover coefficent, c*^ necessary to achieve 0(a) 
improvement and also the optimal critical hopping parameter, k*, eqs. (H5|) . ( H7I) 
over a wide range of coupling constant. 

As a increases we need a significant Csw ^ c*sw^ = 1 for 0{a) improvement. 
We are now seeking a region where a ~ 0.05 — 0.1 fm. Improvement, which is 
presumably represented by an asymptotic series, brings an advantage for smaller 
a say a < 0.1 fm. The two extremes for a are simulations at small a with 'large' 
rrips when there is no continuum extrapolation but a chiral extrapolation, or 
alternatively simulations at 'coarse' a with rups ~ when there is no chiral ex- 
trapolation but a continuum extrapolation. Of course the Schrodinger functional 
does not tell us a; for this conventional HMC simulations are required. Some 
preliminary results indicate that around /3 ~ 5.50 we have a ~ 0.08 fm. 
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Appendix 



A M and AM results 

We collect here in Tables [H El El H El El d E] and [9] the numerical values of M, 
AM as defined in eq. ([27D for the x 2A^, = 8^ x 16 lattices, while in tables [10] 
and [TT] the results for the 12^ x 24 lattices are given. 

The data sets are of size 0(3000) trajectories for the 8^ x 16 lattices and 
0(2000) trajectories for the 12'^ x 24 lattices. An initial thermalisation phase was 
typically of order 300 trajectories. The trajectory length Tchroma was always 1, 
while the number of steps in the trajectory, , , varied for the 8^ x 16 lattices 

d 1 I chroma ' 

from 10 for /? > 6.50 to 12, 12, 15, 18 for /3 = 6.00, 5.50, 5.25, 5.10 respectively. 
This maintained an acceptance rate of > 80%. (This decreased very slightly for 
the larger /5-values.) For the 12^ x 24 lattices rir , = 18, 22 for 13 = 6.00, 5.50 

<~> I / ' chroma ' > ' 

was used to give this acceptance. 

The jackknife errors for the ratios are given uniformly to two significant fig- 
ures, with the overriding requirement that the result must also have a minimum 
of four significant figures. To reduce possible autocorrelations in the data every 
second trajectory was used with a jackknife block size of 10. 



/5 






M 


AM 


5.10 


3.20 


0.11760 


0.007049(2313) 


0.005762(1923) 


5.10 


3.20 


0.11780 


-0.01315(205) 


0.001565(1545) 


5.10 


3.20 


0.11800 


-0.02324(231) 


0.003037(1212) 


5.10 


3.20 


0.11820 


-0.04187(212) 


-0.001353(1782) 


5.10 


3.30 


0.11610 


0.01941(227) 


0.004648(1570) 


5.10 


3.30 


0.11620 


0.001408(2298) 


-0.003942(1737) 


5.10 


3.30 


0.11640 


-0.005654(2058) 


0.001279(1438) 


5.10 


3.30 


0.11660 


-0.02596(166) 


-0.002347(1310) 


5.10 


3.30 


0.11690 


-0.04356(181) 


-0.004137(1550) 


5.10 


3.40 


0.11470 


0.01098(191) 


-0.003299(1305) 


5.10 


3.40 


0.11490 


-0.004606(1516) 


-0.004438(1044) 


5.10 


3.40 


0.11510 


-0.01742(160) 


-0.006135(1442) 


5.10 


3.40 


0.11530 


-0.02432(125) 


-0.005855(748) 


5.10 


3.40 


0.11550 


-0.03424(165) 


-0.004780(1086) 



Table 1: 8^ x 16 results for M and AM for /3 = 5.10. 
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/? 


CgiD 




M 


AM 


5.25 


2.90 


0.12000 


0.02772(322) 


0.007506(2096) 


5.25 


2.90 


0.12015 


0.01468(226) 


0.008527(1944) 


5.25 


2.90 


0.12025 


0.005850(2967) 


0.003412(1559) 


5.25 


2.90 


0.12050 


-0.02097(186) 


0.004020(1167) 


5.25 


2.90 


0.12100 


-0.04947(241) 


0.0008952(14184) 


5.25 


3.00 


0.11860 


0.02041(279) 


-0.0007319(14965) 


5.25 


3.00 


0.11875 


0.0008556(19694) 


0.001173(1115) 


5.25 


3.00 


0.11890 


-0.006210(2160) 


-0.001295(1424) 


5.25 


3.00 


0.11905 


-0.01727(244) 


-0.006479(2808) 


5.25 


3.00 


0.11920 


-0.03280(169) 


-0.002655(1102) 


5.25 


3.10 


0.11700 


0.01973(153) 


-0.001642(991) 


5.25 


3.10 


0.11720 


0.01021(171) 


-0.002551(1054) 


5.25 


3.10 


0.11740 


-0.002194(1506) 


-0.001440(922) 


5.25 


3.10 


0.11760 


-0.01303(133) 


-0.002512(1050) 


5.25 


3.10 


0.11780 


-0.02344(186) 


-0.0008732(12753) 


5.25 


3.20 


0.11580 


0.01019(131) 


-0.005591(815) 


5.25 


3.20 


0.11600 


0.0001673(11516) 


-0.005485(875) 


5.25 


3.20 


0.11620 


-0.008058(1185) 


-0.005259(1342) 


5.25 


3.20 


0.11640 


-0.01905(114) 


-0.003621(1214) 



Table 2: 8^ x 16 results for M and AM for (3 = 5.25. 







K 


M 


AM 


5.50 


2.50 


0.12300 


0.02608(208) 


0.005685(1134) 


5.50 


2.50 


0.12320 


0.01112(215) 


0.003630(1484) 


5.50 


2.50 


0.12335 


-0.001449(2014) 


0.004018(1567) 


5.50 


2.50 


0.12360 


-0.01565(210) 


0.007337(1321) 


5.50 


2.60 


0.12170 


0.006007(1703) 


0.0009700(16220) 


5.50 


2.60 


0.12190 


-0.0001614(18320) 


0.001739(1046) 


5.50 


2.60 


0.12210 


-0.01343(170) 


0.0001628(11051) 


5.50 


2.60 


0.12230 


-0.01959(223) 


0.003397(1508) 


5.50 


2.70 


0.12015 


0.01584(149) 


-0.002008(1139) 


5.50 


2.70 


0.12040 


0.002419(1200) 


-0.001062(798) 


5.50 


2.70 


0.12070 


-0.01264(125) 


-0.001321(1175) 


5.50 


2.70 


0.12090 


-0.01831(150) 


-0.001626(915) 



Table 3: 8^ x 16 results for M and AM for /? = 5.50. 
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/5 




K 


M 


AM 


6.00 


2.10 


0.12430 


0.01841(99) 


0.001623(800) 


6.00 


2.10 


0.12460 


0.006443(1084) 


0.001332(753) 


6.00 


2.10 


0.12495 


-0.004446(970) 


0.0006452(7878) 


6.00 


2.10 


0.12520 


-0.01316(107) 


0.003539(970) 


6.00 


2.20 


0.12330 


0.01135(86) 


-0.0007576(5905) 


6.00 


2.20 


0.12355 


0.002234(706) 


-0.0001747(6084) 


6.00 


2.20 


0.12390 


-0.01050(79) 


-0.0008061(7138) 


6.00 


2.20 


0.12420 


-0.02108(79) 


-0.0008650(6771) 


6.00 


2.30 


0.12190 


0.01996(58) 


-0.002989(439) 


6.00 


2.30 


0.12215 


0.009817(838) 


-0.002765(574) 


6.00 


2.30 


0.12240 


0.0001335(7744) 


-0.003061(672) 


6.00 


2.30 


0.12280 


-0.01430(67) 


-0.003268(549) 


6.00 


2.40 


0.12100 


0.01228(69) 


-0.004705(456) 


6.00 


2.40 


0.12120 


0.003415(610) 


-0.005526(586) 


6.00 


2.40 


0.12140 


-0.004357(723) 


-0.004751(540) 


6.00 


2.40 


0.12160 


-0.01066(73) 


-0.004149(657) 



Table 4: 8^ x 16 results for M and AM for (5 = 6.00. 





^sw 




M 


AM 


6.50 


1.80 


0.12550 


0.01994(59) 


0.001612(472) 


6.50 


1.80 


0.12575 


0.01067(59) 


0.001914(457) 


6.50 


1.80 


0.12600 


0.001513(513) 


0.001973(466) 


6.50 


1.80 


0.12650 


-0.01600(55) 


0.002172(496) 


6.50 


1.90 


0.12440 


0.02139(60) 


0.0004039(4011) 


6.50 


1.90 


0.12470 


0.01068(56) 


0.001113(435) 


6.50 


1.90 


0.12495 


0.001754(539) 


0.0003388(5215) 


6.50 


1.90 


0.12520 


-0.007849(601) 


-0.00003026(52328) 


6.50 


2.00 


0.12360 


0.01255(49) 


-0.002074(450) 


6.50 


2.00 


0.12390 


0.001931(525) 


-0.001253(358) 


6.50 


2.00 


0.12410 


-0.006006(505) 


-0.002711(510) 


6.50 


2.00 


0.12110 


-0.01635(19) 


-0.001291(153) 



Table 5: 8^ x 16 results for M and AM for /? = 6.50. 
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P 






M 


A M 


7 9(1 


1 ZLO 
1 .^U 


n 1 9790 


0.02534(46) 




u.uuuouo^^oo ( j 


7 90 


1 40 
1 .^u 


n 1 97Q7 


-0.0007597(4109) 




00'^09Q('4'^0'l 
u .uutjuzy I ^ou I 


7 90 


1 40 
± .^u 




r\ r\-\ '-J -\ / A C>\ 

-0.01713(48) 




00^1 1 

U.UUtJllOlO 1 Z I 


7 90 


1 40 


19Q9n 


-0.03970(53) 




007053l'563'l 

U.UU 1 U^JOl ^JUO ) 


7 90 


1 fiO 

i .uu 


1 Qi^no 

U. iZOUU 


0.03839(43) 






7 90 


1 fiO 
i .uu 


1 9P;70 
U. iZO ( u 


0.01500(38) 




001 p;'?4l''?8Q'\ 
u.uuioo^^^ooy j 


7 90 


1 fiO 
1 .uu 


1 9(S1 ^ 
u. 1 zu 1 tj 


0.0003391(4484) 




001 SS'^l'^^l 1 
u.uuiooo^oii J 


7 90 
( . zu 


1 fiO 
1 .uu 


1 9(S(S0 

U. IZUUU 


-0.01525(36) 




001 i^k^aXx 

U.UUlOOOi^O^O j 


7 90 
( . zu 


1 (SO 
1 .uu 


1 9790 
u. xz 1 zu 


-0.03608(38) 




001 PiAACMM^ 

U.UUXU^^lOUl 1 


7 90 
( .zu 


1 80 
i .oU 


1 9970 
u. izz / u 


0.05607(29) 




001 7S(S('9'?Q'\ 
-u.uui / ou^^zoy j 


7 90 
( .ZU 


1 80 
i .oU 


1 9'?S0 

U. IZOOU 


0.01959(32) 




001 '^i^'?i'9(S0'i 

-U.UUlOOd^^ZUU J 


7.20 


1.80 


12438 


-0.00008070(34186) 


-0 002103f288) 


7.20 


1.80 


0.12500 


-0.02136(34) 




-0.001939(300) 


7.20 


1.80 


0.12590 


-0.05319(35) 




-0.001455(388) 


7.20 


2.00 


0.12150 


0.03819(31) 




-0.005604(315) 


7.20 


2.00 


0.12210 


0.01736(38) 




-0.005245(340) 


7.20 


2.00 


0.12264 


-0.0002027(3196) 




-0.005262(470) 


7.20 


2.00 


0.12290 


-0.008518(356) 




-0.005990(375) 


7.20 


2.00 


0.12360 


-0.03421(34) 




-0.006188(311) 


Table 6: 8^ x 16 results for M and AM for (5 = 7.20. 






K 


M 


AM 


8.00 


1.40 


0.12570 


0.02742(25) 


0.003117(207) 


8.00 


1.40 


0.12630 


0.007469(239) 


0.002932(272) 


8.00 


1.40 


0.12651 


0.0001971(2329) 


0.002716(221) 


8.00 


1.40 


0.12680 


-0.009671(223) 


0.003270(247) 


8.00 


1.40 


0.12730 


-0.02596(28) 


0.003221(256) 


8.00 


1.60 


0.12430 


0.02266(23) 


-0.0004972(2019) 


8.00 


1.60 


0.12480 


0.005679(245) 


-0.0008718(2676) 


8.00 


1.60 


0.12498 


0.0002484(2335) 


-0.0008608(2491) 


8.00 


1.60 


0.12520 


-0.007169(242) 


-0.0006004(2378) 


8.00 


1.60 


0.12570 


-0.02410(25) 


-0.001201(239) 


8.00 


1.80 


0.12240 


0.03501(24) 


-0.003785(264) 


8.00 


1.80 


0.12290 


0.01858(26) 


-0.003763(179) 


8.00 


1.80 


0.12344 


0.0005959(2472) 


-0.004154(247) 


8.00 


1.80 


0.12350 


-0.002196(264) 


-0.005071(223) 


8.00 


1.80 


0.12400 


-0.01861(27) 


-0.004060(270) 



Table 7: 8^ x 16 results for M and AM for /3 = 8.00. 
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/3 




K 


M 


AM 


10.00 


1.20 


0.12570 


0.01641(22) 


0.003409(179) 


10.00 


1.20 


0.12619 


-0.0001306(1605) 


0.003338(182) 


10.00 


1.20 


0.12630 


-0.003541(173) 


0.003321(217) 


10.00 


1.20 


0.12690 


-0.02350(20) 


0.003296(198) 


10.00 


1.20 


0.12750 


-0.04340(17) 


0.003247(206) 


10.00 


1.40 


0.12410 


0.03094(21) 


0.0003442(1695) 


10.00 


1.40 


0.12470 


0.01351(29) 


0.0006860(2239) 


10.00 


1.40 


0.12507 


0.0004563(4134) 


0.001032(171) 


10.00 


1.40 


0.12530 


-0.008549(319) 


0.0005683(2022) 


10.00 


1.40 


0.12590 


-0.02794(27) 


0.001172(222) 


10.00 


1.60 


0.12270 


0.03342(46) 


-0.004086(267) 


10.00 


1.60 


0.12320 


0.02152(16) 


-0.003744(145) 


10.00 


1.60 


0.12382 


0.001171(165) 


-0.003759(157) 


10.00 


1.60 


0.12390 


-0.002455(294) 


-0.003601(186) 


10.00 


1.60 


0.12450 


-0.02161(19) 


-0.004090(163) 



Table 8: 8^ x 16 results for M and AM for P = 10.00. 









M 


AM 


14.00 


1.10 


0.12500 


0.01723(8) 


0.002646(110) 


14.00 


1.10 


0.12530 


0.007452(89) 


0.002787(103) 


14.00 


1.10 


0.12560 


-0.002273(87) 


0.002941(114) 


14.00 


1.10 


0.12590 


-0.01196(9) 


0.002676(100) 


14.00 


1.10 


0.12620 


-0.02194(9) 


0.002684(113) 


14.00 


1.20 


0.12420 


0.03218(36) 


0.001696(167) 


14.00 


1.20 


0.12470 


0.01423(16) 


0.001044(98) 


14.00 


1.20 


0.12530 


-0.002225(329) 


0.002191(174) 


14.00 


1.20 


0.12580 


-0.01786(46) 


0.002320(199) 


14.00 


1.30 


0.12380 


0.02900(34) 


-0.001514(170) 


14.00 


1.30 


0.12430 


0.01132(40) 


-0.0004894(1641) 


14.00 


1.30 


0.12480 


-0.005572(301) 


-0.0007392(1529) 


14.00 


1.30 


0.12530 


-0.02027(107) 


-0.0009807(2838) 



Table 9: 8^ x 16 results for M and AM for P = 14.00. 
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M 


AM 


5.50 


2.50 


0.12300 


0.02540(221) 


0.0006153(17988) 


5.50 


2.50 


0.12320 


0.004367(3139) 


0.001051(1568) 


5.50 


2.50 


0.12335 


-0.002279(2162) 


0.001425(1272) 


5.50 


2.50 


0.12360 


-0.01981(151) 


0.002050(1237) 


5.50 


2.60 


0.12170 


0.007744(2026) 


-0.0009438(8558) 


5.50 


2.60 


0.12190 


-0.002810(1805) 


0.0002407(11134) 


5.50 


2.60 


0.12210 


-0.01117(179) 


-0.0008471(12790) 


5.50 


2.60 


0.12230 


-0.02560(168) 


-0.0007416(19293) 


5.50 


2.70 


0.12015 


0.01289(175) 


-0.0008967(10745) 


5.50 


2.70 


0.12040 


0.001170(2838) 


-0.0003062(16169) 


5.50 


2.70 


0.12070 


-0.01224(140) 


-0.001362(852) 


5.50 


2.70 


0.12090 


-0.02138(153) 


0.0002677(9645) 



Table 10: 12^ x 24 results for M and AM for /3 = 5.50. 





Csw 


K 


M 


AM 


6.00 


2.10 


0.12430 


0.01957(74) 


0.0003629(5316) 


6.00 


2.10 


0.12460 


0.007496(680) 


0.0006202(5838) 


6.00 


2.10 


0.12495 


-0.001642(1038) 


0.001463(1070) 


6.00 


2.10 


0.12520 


-0.01123(113) 


0.0005411(5241) 


6.00 


2.20 


0.12330 


0.01228(67) 


-0.0008308(5383) 


6.00 


2.20 


0.12355 


0.002046(917) 


-0.0008953(4855) 


6.00 


2.20 


0.12390 


-0.01153(83) 


0.0005139(5375) 


6.00 


2.20 


0.12420 


-0.02019(76) 


-0.0003129(6525) 


6.00 


2.30 


0.12190 


0.02111(49) 


-0.001234(455) 


6.00 


2.30 


0.12215 


0.01067(68) 


-0.001233(833) 


6.00 


2.30 


0.12240 


0.002555(557) 


-0.0008735(5407) 


6.00 


2.30 


0.12280 


-0.01306(64) 


-0.0001565(5009) 


6.00 


2.40 


0.12100 


0.01273(49) 


-0.001217(461) 


6.00 


2.40 


0.12120 


0.005458(635) 


-0.002194(415) 


6.00 


2.40 


0.12140 


-0.003718(533) 


-0.002257(514) 


6.00 


2.40 


0.12160 


-0.009398(475) 


-0.001493(486) 



Table 11: 12^ x 24 results for M and AM for (3 = 6.00. 
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